Color transparency and suppression of high-px hadrons in nuclear collisions 
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The production length Ip of a leading (large Zh) hadron produced in hadronization of a highly 
virtual high-pT parton is short because of the very intensive vacuum gluon radiation and dissipation 
of energy at the early stage of process. Therefore, the main part of nuclear suppression of high-pT 
hadrons produced in heavy ion collisions is related to the survival probability of a colorless dipole 
propagating through a dense medium. This is subject to color transparency, which leads to a steep 
rise with pr of the nuclear ratio RAAipr), in good agreement with the recent data from the ALICE 
experiment at LHC, CERN. No adjustment, except the medium density is made, and the transport 
coefficient is found to be qo = 0.8 GeV'^/ fm. This is close to the value extracted from the analysis 
of RHIC data for J/'i' suppression, but is an order of magnitude smaller than the value found from 
jet quenching data within the energy loss scenario. Although the present calculations have a status 
of a postdiction, the mechanism and all formulas have been published, and are applied here with 
no modification, except the kinematics. At the same time, pT-dependence of Raa at the energy of 
RHIC is rather fiat due to the suppression factor steeply falling with rising xt, related to the energy 
conservation constraints. This factor is irrelevant to the LHC data, since xt is much smaller. 

PACS numbers: 24.85.+p, 25.75.-q, 25.75.Bh, 25.75.Cj 



I. INTRODUCTION 

One of the first results of the heavy ion program at 
LHC is the observation by the ALICE experiment [1] of 
a strong nuclear suppression of high-pr charged hadrons. 
These data expose novel features compared with similar 
measurements at RHIC First, the nuclear suppres- 

sion factor Raa reaches significantly smaller values. This 
is not a surprise, since at the LHC energies hadrons orig- 
inate mainly from hadronization of gluons, which have 
a larger color charge than quarks dominating at RHIC. 
Correspondingly, gluons dissipate energy with a higher 
rate. Second, Raa{pt) steeply rises with pT, while it 
exposes a rather flat pT-dependence in RHIC data. The 
latter is affected by the restrictions imposed by energy 
conservation [j . It was predicted Q that the production 
rate for hadrons and direct photons is suppressed in pA 
and AA collisions by the deficit of energy not only at 
forward rapidities, but also at large xt — '2pt/\/s- Here 
we concentrate on the interpretation of LHC data, which 
are free of these complications, since the values of xt are 
very small. 

The popular model explaining the observed suppres- 
sion of Raa at high-py relates it to the induced radiation 
energy loss by a parton propagating through the medium, 
which was created in the nuclear collision (e.g. see in 0). 
This energy-loss scenario is based on the unjustified as- 
sumption that hadronization of the parton lasts longer 
than the time of propagation through the medium, and 
that the detected hadron is always produced outside the 
medium. However, because of the steeply falling px de- 
pendence of the cross section, most of high-pT hadrons 
carry a large fraction of the jet momentum, and en- 
ergy conservation constraints the production length Ip for 



such leading hadrons. It is expected to be rather short 
even within the simple string model 0, and should 
be much shorter in the case of intensive vacuum gluon 
radiation by a high-p^ parton. 

One should clearly distinct between the production 
time scales for a colorless dipole (pre-hadron) and the 
final hadron. The former signals on color neutralization, 
which stops the intensive energy loss caused by vacuum 
radiation following the hard process, while the latter is a 
much longer time taken by the dipole to gain the certain 
hadronic mass, i.e. to develop the hadron wave func- 
tion. While the former contracts oc {1 — Zh) at large 
fractional momentum of the hadron, the latter keeps 
rising cx z^- These two time scales are frequently mixed 
up. The shortness of the production lengths at large Zh is 
dictated by energy conservation. Indeed, a parton orig- 
inated from a hard reaction intensively radiates losing 
energy, and this should cannot last long, otherwise the 
parton energy will drop below the energy of the detected 
hadron. 

One should also distinct between the mean hadroniza- 
tion time of a jet, whose energy is shared between many 
hadrons, and specific events containing a leading hadron 
with z/i — )• 1. Production of such a hadron in a jet is 
a small probability fiuctuation, usually associated with 
large rapidity gap events. The space-time development 
of such an unusual jet is different from the usual averaged 
jet. 

The controversy between the models with short and 
long production times has been under debate the last two 
decades (see in 0, Q), but no proof of a long time scale 
has been proposed so far, to the best of our knowledge. 
Of course, an experimental verification would be most 
convincing. Data on high-pr hadron production in heavy 



ion collisions provides a rather poor test of the models. 
Too many uncertainties are involved, the medium proper- 
ties are unknown, and their variation in space and time is 
based on simplified, even ad hoc models. The fractional 
energy is not known, but enters the convolution of 
the initial parton distribution, hard cross section and the 
fragmentation function. The important contribution of 
initial state effects (cold nuclear matter) can be only cal- 
culated within models. 

Probably the best way to study the space-time develop- 
ment of hadronization is inclusive hadron production in 
DIS on nuclei at large Bjorken x. In this case the medium 
density and its spacial distribution is well known. The 
fractional energy Zh of the hadron is directly measured. 
A good model should predict the nuclear modification 
factor with no fitting. Indeed, such a prediction was pro- 
vided in [ll| within the hadronization model with a fi- 
nite production length. Later the first data from the 
HERMES experiment [12,] confirmed well this prediction. 
Also HERMES results on broadening of transverse mo- 
mentum, sensitive to the production length were ex- 
plained well in [lo| . On the other hand, the compre- 
hensive study of nuclear effects within the pure energy 
loss scenario performed recently in |13j led to striking 
disagreement with data for leading hadrons. 

Here we rely on the model [14] for the production time 
distribution of leading hadrons in a jet, produced at the 
mid rapidity. In this case the initial parton energy and 
virtuality are equal. 



E = Q = kT = —, 

Zh 



(1) 



where kr and px are the transverse momenta of the par- 
ton initiating the jet and of the detected hadron, respec- 
tively. An example of the Ip distribution at z^ = 0.7 and 
different quark jet energies is shown in Fig. [1] One can 



0.3 



0.2 



Q 

T3 



0.1 





1 1 1 1 ■ 1 ■ 1 1 

z,=o.7 : 




E=6 GeV ; 


E=100GeV^ 





10 



Ip (fm) 



FIG. 1: The pre-hadron production length distribution 
dD{z) /dip (in arbitrary units) for for a quark jet with en- 
ergies E = kr ^ 6, 10, 20, 100 GeV (from top to bottom) 
and Zh = 0.7. 

see that the Zp-distribution narrows with energy, but is 



leveling off at high energies. This happens due to com- 
pensation of several effects [l3] , acting in opposite direc- 
tions. The Lorentz factor makes Ip longer with energy, 
while the increasing virtuality gives rise to a more inten- 
sive gluon radiation and energy loss in vacuum, leading 
to a shorter Ip. Moreover, the Sudakov suppression, es- 
sential at large Zh, also shortens Ip. 

The mean value (Ip) for quark jets, weighted with the 
distribution dD/dlp, is depicted by solid curves in Fig.[2l 
as function of energy, for Zh = 0.5, 0.7, 0.9. Indeed, {Ip) 
saturates at high jet energies E = kx. 

Notice that in DIS the mean production length at 
fixed rises linearly as function of energy. This case 
is quite different from high-py, where the virtuality in- 
creases with energy. 

For gluon jets the energy loss is larger and the Su- 
dakov suppression stronger due to the Casimir factor. 
This leads to a shorter production length. The mean 
value of Ip for gluon jets as function of energy is depicted 
in Fig. [2]by dashed curves, for several values of Zh. We 
see that at high energies {Ip) < Ifm. Thus, one can say 
that the pre-hadron is created almost instantaneously, 
because its production time is shorter than the expected 
time of medium creation, to ~ 1 fm- 

An interesting possibility of a very dense medium was 
considered in [5|, [l^ . If the medium were so dense that 
the mean free path of the produced dipole was vanish- 
ingly small, the nuclear suppression factor Raa would be 
proportional to {Ip), and could be predicted in a parame- 
ter free way. Although such a possibility does not contra- 
dict RHIC data on jet quenching 0, [ll] , the recent study 
of J/^* suppression in heavy ion collisions [13, [3 found 
that the density of the medium is rather low. The max- 
imal value of transport coefficient reached at the time 
scale of the medium creation <o ~ 1 fni was found to 
be quite low, go ~ 0.2GeV^/fm, an order of magnitude 
smaller than follows from hadron suppression data in- 
terpreted within the energy loss formalism p^ . Thus, 




FIG. 2: The mean production length as function of energy for 
quark (solid curves) and gluon (dashed curves) jets. In both 
cases the curves are calculated at Zh = 0.5, 0.7, 0.9 (from top 
to bottom). 
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the scenario of a very dense medium is not supported 
by data, and one should study propagation of colorless 
dipole through a moderately opaque medium. 

Here we employ the same description of dipoles in a 
medium, as was used for J/^ suppression [17|, and apply 
it to high-px processes. The main difference is the time 
scale of formation of the final hadron: while it is very 
short for a heavy and slow J /'^ [13, [13; the formation 
of the wave function of a light qq dipole moving with 
a high momentum is rather long. Correcting for this, 
we apply the same description to the suppression nuclear 
factor Raa{pt) for hadrons produced in central lead-lead 
collisions at LHC, and achieve good agreement with the 
recent data from the ALICE experiment [1]. Moreover, 
we found a value of the transport coefficient similar to 
what was found from J/\l/ data at RHIC. 



II. ATTENUATION OF A SMALL SIZE DIPOLE 
IN A MEDIUM 

A parton, which experienced a hard scattering with 
transverse momentum fcy, shakes off its color field in the 
form of a forward cone of gluon radiation, and starts 
propagating along a new direction, lacking the soft part 
of its field up to transverse frequencies < fc^. The parton 
starts regenerating its color field, with transverse size r(/) 
expanding as function of path length starting from the 
initial small size tq '-^ 1/fcr at / = 0. At some distance Ip 
a colorless dipole, pre-hadron, is produced with a size of 
the order of the transverse size r{lp) of the regenerated 
field. We describe the regeneration grocess and gluon 
radiation within the dipole approach [19, 20]. 

A small size dipole is expanding so fast that its initial 
size is quickly forgotten. Indeed, the speed of expansion 
of a dipole correlates with its size: the smaller the dipole 
is, the faster it is evolving. This is controlled by the 
uncertainty principle, q ~ l/r. 
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where E — pT '\s the dipole energy in the cm. of the 
collision; vt and q ^ l/r are the transverse velocity and 
momentum of the quark relative to the dipole momentum 
direction. The solution of this equation reads [13, [3 j 
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where vq l/pr is the initial dipole size, neglected in 
what follows. 

Propagation of a dipole over path length i in a medium 
is characterized by a survival probability 



S{L) = exp 



dl a[r{l)] p{l) 



(4) 



where the dipole cross section a(r) times the medium 
density p is the attenuation rate of the dipole. 



The dipole cross section for small dipoles is a{r) = 
C . The factor C, for dipole-proton interactions, is 
known from DIS data. Its value for a hot medium is un- 
known, as well as the medium properties. It is convenient 
to express it in terms of the so called transport coefficient, 
which is the broadening of a parton in the medium over 
the path length 1 fm. Indeed, the same factor C controls 
both the dipole cross section and broadening of a quark 
propagating through the medium ^ll[2^. So, the factor 
C in the dipole cross section is related to the transport 
coefficient q [2^ , which is the in- medium broadening per 
unit of length. 



2p 
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Then, using ([3]) (with t = I) one can represent the sur- 
vival probability of the dipole in the medium, Eq. (j4|. 

as 



S{L) = exp 



dlq{l)r^{l) 



exp 



1^ 

-— [ dllqil) 
PT J 



(6) 



Now we are in a position to calculate the nuclear atten- 
uation factor for a high pT hadron produced in heavy ion 
collision. For central (6 = 0) collisions of two identical 
nuclei one should integrate over the impact parameter t 
of the hard collision, with a weight factor T^(t), where 
Ta(t) = J_^dz pA{f, z) is the nuclear thickness func- 
tion, the integral of nuclear density along the collision 
direction: 

oo 

Jd^TTlir) Raa{t,Pt) 
RAAib = 0,pT) = ^ . (7) 



The factor Raa{t,Pt) is the nuclear suppression factor 
corresponding to production of a high-fc^ parton at im- 
pact parameter r, propagating then over a path length 
(Ip), radiating gluons and losing energy, and eventually 
producing a colorless dipole pre-hadron with transverse 
momentum px = kx / zu, which propagates through the 
nucleus evolving its size according to Eq. ([3]). We rely on 
the above evaluation of {Ip) in vacuum, since the medium- 
induced energy loss is much smaller than the vacuum one. 
Besides, induced energy loss can make Ip only shorter, 
which will not affect the further calculations. This sup- 
pression factor has the form [13, [11] , 
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maxllpJo}. Here = 1 fm is the adjusted go = 1.6GeV^/fm. Since at large > 100 GeV 



where l^a 

time scale of creation of the medium resulted from gluon 
radiation at mid rapidities in heavy ion collisions. Since 
the production length for a gluon jet is short, (Ip) < Iq, 
its actual value is not important. 

The medium density is time dependent, and is assumed 
to dilute as p{t) = poto/t due to the longitudinal ex- 
pansion of the produced medium. Correspondingly, the 
transport coefficient depends on impact parameter and 
time (path length) as [24| . 



q{l, b, f) 



qolo npart{b,T) 



I "part (0,0) 



(9) 



where npart(b,f) is the number of participants; cor- 
responds to the maximum medium density produced at 
impact parameter r = in central collisions {b — 0) at 
the time t = to — lo after the collision. In what fol- 
lows we treat the transport coefficient go as an adjusted 
parameter. 



III. RESULTS VS DATA 

The results of the calculation with Eqs. ©-([5]) with 
qa = 0.8 GeV^/ fm are shown by the solid curve in Fig. [31 
in comparison with ALICE data T|. Except for go, no 



valence quarks with larger Ip should dominate, the rise 
of Raa with pt should slightly slow down deviating from 
the solid curve towards the dashed one. 

We do not attempt here to describe Raa at pt < 
7 GeV, since the dynamics becomes much more compli- 
cated. First, the production length becomes several time 
longer, as is depicted in Fig. [21 Second, one should take 
into account the Cronin effect in A A collisions, which is 
poorly known because of the large fraction of baryons in 
detected charged hadrons. 

For the same reason, application of this mechanism 
to high-pT hadron production at RHIC goes beyond the 
scope of this paper, since it involves a more complicated 
and model dependent dynamics. Indeed, it was demon- 
strated in that for particle production at forward ra- 
pidities, Xp ^0.1, energy conservation becomes an issue. 
It causes additional nuclear suppression, steeply increas- 
ing with xp. The same, of course, should happen at 
large xt = 2pT / y/s. It was demonstrated in fS] that this 
mechanism leads to suppression of high-p^^ hadrons and 
direct photons in pA and A A collisions. Therefore, be- 
sides other reasons (different kinematics, valence quark 
dominance, etc.), this mechanism makes the pT depen- 
dence of Raa{pt) significantly flatter at RHIC than at 
LHC. Detailed calculations for the RHIC energy domain 
will be published elsewhere. 




Pt (GeV) 



FIG. 3: ALICE data for central, 0-5%, lead-lead collisions at 
2.76 GeV [ij. The solid curve corresponds to Eqs. 



([7|) calculated with Ip < Iq = lira and go 
The dashed curve, 
sensitivity to Ip. 



calculated with L = 2fm, 



O.SGeVVfm. 
demonstrates 



further adjustment was done, and agreement with data at 
large Pt > 7 GeV is pretty good. Moreover, the transport 
coefficient turns out to be of the same order as was found 



17| from data on J/^E" production. 



The dashed curve in Fig. [31 demonstrates the sensitivity 

2 fm and accordingly 



to Ip. It is calculated at {Ip 



IV. SUMMARY 

Although these calculations have a status of post- 
diction, the mechanism and all formulas have been al- 
ready published, and we apply them here with no spe- 
cific modification, except for the kinematics. The dy- 
namics of nuclear suppression of high-py hadrons pro- 
duced in central lead-lead collisions at = 2.76 GeV is 
based on the shortness of the production length of a pre- 
hadron, and its development and propagation through 
a dense medium. We performed calculations within the 
same scheme as was used for the analysis of J/^' 
production data. Moreover, we arrived at a value of the 
transport coefficient, characterizing the properties of the 
medium, which is pretty close to the value extracted from 
J/4' data, and is an order of magnitude smaller than 
what has been obtained from analyses of jet quenching 
RHIC data based on the energy loss scenario. 

These calculations can be improved by replacing the 
simplified description of the dipole evolution, Eq. ([5t . 
with the rigorous quantum-mechanical approach [25| 
based on the path integral technique. This description 
can be also applied to RHIC data on jet quenching, but 
one should introduce a model dependent suppression fac- 
tor, which is related to the constraints on nuclear parton 
distributions imposed by energy conservation [4]. These 
further developments of the present approach will be pub- 
lished elsewhere. 



4 



Acknowledgments 

We are grateful to Jan Nemchik for pointing out a 
dimension related typo in the code. We are also thankful 
to Martin Poghosyan who provided us with the tables 



of ALICE data. This work was supported in part by 
Fondecyt (Chile) grants 1090236, 1090291 and 1100287, 
by DFG (Germany) grant PI182/3-1, and by Conicyt- 
DFG grant No. 084-2009. 



[1] K. Aamodt et al. [ALICE Collaboration], 

larXiv:1012.1004l [nucl-ex]. 
[2] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. 

Lett. 91, 072301 (2003). 
[3] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 

91, 072304 (2003). 
[4] B.Z. Kopeliovich, J. Nemchik, I.K. Potashnikova, 

M.B. Johnson,and I. Schmidt, Nucl. Phys. B146, 171 

(2005). 

[5] B. Z. Kopeliovich and J. Nemchik, arXiv:1009.1162l [hep- 
ph]. 

[6] S. Wicks, W. Horowitz, M. Djordjevic and M. Gyulassy, 

Nucl. Phys. A 784, 426 (2007). 
[7] B.Z. Kopeliovich and F. Niedermayer, Sov. J. Nucl. Phys. 

42 (1985) 504; Yad. Fiz. 42 (1985) 797. 
[8] A. Bialas and M. Gyulassy, Nucl. Phys. B 291, 793 

(1987). 

[9] B. Z. Kopeliovich, J. Nemchik, E. Predazzi and 
A. Hayashigaki, Nucl. Phys. A 740, 211 (2004). 

[10] S. Domdey, D. Grunewald, B. Z. Kopeliovich and 
H. J. Pirner, Nucl. Phys. A 825, 200 (2009). 

[11] B.Z. Kopeliovich, J. Nemchik and E. Predazzi, Proceed- 
ings of the workshop on Future Physics at HERA, ed. 
by G. Ingelman, A. De Roeck and R. Klanner, DESY 
1995/1996, V. 2, 1038; arXiv: nucbth/9607036 

[12] A. Ahapetian et al. [HERMES Collaboration], Eur. 
Phys. J. C 20, 479 (2001). 

[13] W. t. Deng and X. N. Wang, Phys. Rev. C 81, 024902 



(2010). 

[14] B. Z. Kopeliovich, H. J. Pirner, I. K. Potashnikova and 

I. Schmidt, Phys. Lett. B 662, 117 (2008). 
[15] N. Armesto et al, J. Phys. G 35, 054001 (2008); 

B. Z. Kopeliovich, I. K. Potashnikova and L Schmidt, 

arXiv:0707.4302 [nucl-th]. 
[16] A. Adare et al. [PHENIX Collaboration], Phys. Rev. C 

77, 064907 (2008). 
[17] B. Z. Kopeliovich, I. K. Potashnikova and I. Schmidt, 

Phys. Rev. C 82, 024901 (2010). 
[18] B. Z. Kopeliovich, arXiv:1007.4513 [hep-ph]; to appear 

in Nucl. Phys. A. 
[19] B. Z. Kopeliovich, A. V. Tarasov and A. Schafer, Phys. 

Rev. C 59, 1609 (1999). 
[20] B.Z. Kopeliovich, A. Schafer and A.V. Tarasov, Phys. 

Rev. D62 (2000) 054022. 
[21] M. B. Johnson, B. Z. Kopeliovich and A. V. Tarasov, 

Phys. Rev. C 63, 035203 (2001). 
[22] J. Dolejsi, J. Hiifner and B. Z. Kopeliovich, Phys. Lett. 

B 312, 235 (1993). 
[23] R. Baier, Y. L. Dokshitzer, S. Peigne and D. Schiff, Phys. 

Lett. B 345, 277 (1995). 
[24] X. F. Chen, C. Greiner, E. Wang, X. N. Wang and Z. Xu, 

Phys. Rev. C 81, 064908 (2010) 
[25] B. Z. Kopeliovich and B. G. Zakharov, Phys. Rev. D 44, 

3466 (1991). 



5 



